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NOTATION 
Cross sectional area of towing tank 
One-half of length of prolate spheroid 
Beam of ship 
Radius of maximum cross section of spheroid 


Block coefficient 


Frictional resistance coefficient 


Frictional Resistance Coefficient 
Total resistance coefficient 
Wave resistance coefficient 
Prismatic coefficient 


Hull fineness parameter 


Water depth Froude number 
Ship length Froude number 


Acceleration of gravity 
Water depth 
Potential jump due to blockage 


Partial form factor 


Body length; length between perpendiculars 


Length of waterline 


Added mass of submerged body 


Added mass coefficient 


Ww 


Radius in cylindrical coordinates 


R, Frictional resistance 
Rp Total resistance 


Wave resistance 


RY Reynolds number 

Ro Radius of tunnel wall 

Sy Wetted surface 

T Draft of ship 

U Uniform incoming stream velocity at upstream 
infinity 

u Mean speed due to blockage 

uU Mean speed on the body in unbounded water 

W Width of towing tank 

SVs Right-handed rectangular coordinates 

Au Speed increment 

Au Mean-speed increment averaged over body 

v Kinematic viscosity of water 

fo) Density of water 

T= (T) T5573) Tangential unit vector 

co) Total velocity potential 

oe Total velocity potential in absence of 


tank (or tunnel) walls 
0) Perturbation velocity potential 


0) Perturbation velocity potential in absence of 
tank (or tunnel) wall 


vi 


yy 


V 


Displaced volume; volume 


Gradient operator 


= 
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ABSTRACT 


It has recently been shown that a jump in velocity potential 
exists between infinite upstream and downstream directions when a 
body translates unformly along a channel of finite cross section 
such as a towing tank or wind tunnel. In this report a new block- 
age correction formula for body speed is proposed. The speed 
correction formula due to blockage is obtained by dividing com- 
puted potential jump by body length, assuming that the body is 
slender or flat in the direction of motion. The potential jump 
is expressed explicitly in terms of the effective volume, i.e., 
the sum of the displaced volume and added mass/density of the 
submerged body, and the depth Froude number, if a free surface 
is present. As a test of the present speed correction formula, 
two cases are considered: (1) the Wigley parabolic ship model, 
tested in both a small and a large towing tank, (2) a body of 
revolution (prolate spheroid) tested in a circular wind tunnel. 
In each case the mean-speed increment averaged over the entire 
body surface is computed by a three-dimensional, finite-element 
method applicable to free-surface flow problems. These are shown 
to be in good agreement with those obtained by the approximate 
speed correction formula. At high values of Froude numbers, the 
main difference in the total resistance coefficients measured in 
the two towing tanks by Tamura is due primarily to difference in 
model wave resistance computed for the two tanks by a full- 
fledged, three-dimensional, finite-element method. Results are 
also compared to those obtained by using the speed correction 
formula of Lock and Johansen. The present formula renders a 
better approximation than that of Lock and Johansen when the 
cross sectional area of a flow tunnel is not much larger than 
the maximum cross section area of the body. 
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INTRODUCTION 


Many authors have investigated blockage effect and proposed approxi- 


mate blockage formulas to account for towing tank or wind tunnel bound- 


aries. 


1=9* 


The first approximation concerning towing tank blockage effects 


date back more than four decades. Owing to the difficulty encountered in 


computing flow separation, wake flow, free-surface effects, etc., the exact 


*A complete listing of references is given on page 21. 
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magnitude of the blockage effect on fluid force acting on a body is too 
complicated to analyze by purely theoretical means. However, these diffi- 
culties did not stop engineers from attempting to make simple engineering 
approximations of the blockage problem. For engineering purposes, compu- 
tation of a mean-speed increment on a body due to blockage effects has been 
the main focus of interest in order to make a blockage correction to 
frictional drag. In the computation, the incremented change in frictional 
drag due to blockage is determined directly from the computed incremental 
increase of mean speed over the body surface caused by flow blockage. 

Two basic inviscid flow-theory approaches have been previously em- 
ployed. The first approach is based on the so-called one-dimensional, 
mean-flow theory, using the Kreitner equation, which was first obtained by 
eretener from the Bernoulli and the mass continuity equations under the 
assumption that velocity is uniform in each cross sectional plane. To 
name a few, epee and Kim? used this approach. The second approach is 
based on successive reflection of images in the walls of a rectangular tank 
or simpler axisymmetric singularities in case of axisymmetric flows. In 
this approach, the velocity potential of the flow inside a specified tank 
boundary can be computed exactly in principle; usually, the potential is 
represented by a series expansion, and only the first few terms are 


>" and Landweber and Naleayemae have used the 


computed. Oniigeraa,© Tamura, 
latter approach. 

In all, there exist about a dozen formulas proposed for blockage 
corrections and each is somewhat different from the other. Some formulas 
introduce empirical correction faeroren whereas others claim to be based 
on analytical derivations. Some formulas are proposed to be used only for 
frictional resistance corrections, whereas other formulas are used for 
total resistance corrections. An extensive review of the subject has been 
made by Gross and Watanabenc 

In the present preliminary study, skepticism is exercised about 
proposals in speed correction formulas that can be used to correct the 


total resistance which include the wave resistance in water of finite depth 


with sidewalls. Herein is proposed a new speed correction formula to be 
used only for frictional resistance. The wave resistance which has been 
computed for any towing tank and/or model conditions by using the localized 
finite-element method previously developed by the em gnerey It seems 

to be impossible to make blockage corrections to total resistance by using 
only a single-speed correction formula, even though such formulas have 
been proposed in the past. 

The approach used to derive a blockage correction formula herein is 
different from the two inviscid flow-theory approaches described pre- 
viously. Derivation of a mean-speed correction formula in this report is 
based on the potential jump occurring in the three-dimensional flow in a 
towing tank or wind tunnel. To test the new speed correction formula, 
numerical computations for a full-fledged, three-dimensional wave resist- 
ance problem were made. The numerical mean-speed increment on a spheroid 
was computed exactly for a circular wind tunnel and compared with the 
results obtained by the new formula; results obtained agree reasonably 


well with exact numerical results. 


BLOCKAGE CORRECTION 
EXACT MEAN-SPEED INCREMENT 
Steady uniform flow past a ship fixed in a channel has been con- 
sidered; see sae? The coordinate system is right handed and rectangular. 
Under the usual assumptions, steady uniform flow may be described by a 


total velocity potential ® defined by 
® (x,y,z) = Ux + > (x,y,z) (1) 


where » is the perturbation-velocity potential in a channel of finite 


cross section. Similarly the total velocity potential 


e. (x,y,z) = Ux + Oe (x,y,z) (2) 


is defined to describe the flow about the same body in an unbounded fluid, 
i.e., in the absence of channel boundaries. The fluid speed on a body 
surface in general increases due to the blockage effect when compared with 
that of unbounded fluid. However, the speed increment on the body surface 
is not uniform over the entire surface. For example, the forward stag- 
nation point of an axisymmetric body remains the same whether in an un- 
bounded fluid or in a wind tunnel of circular cross section. Nevertheless, 
a mean speed correction has been traditionally employed for the blockage 
correction mainly due to its simplicity. To describe a mean-speed: incre- 
ment, speed increment due to blockage locally on the body surface is 


defined as 


Ay 


Au 


V (aon) O 


(3) 


ay 


V (6-5) ° 


> 
where T = (T, 57 T3) is a unit tangential vector on the body surface; T 


Die 1 
is the component along the x-axis, i.e., the longitudinal direction, and 
To and T3 are, respectively, the normal and tangential components in the 
cross sectional plane of the body. Then the "exact'' mean speed increment 


averaged over the entire submerged body surface is given by 


Agee - {J V (6-9,) * T ds (4) 
S 
(e) 


oO 


where 55 is the wetted surface area, and Tt is specified. One natural way 
of specifying Tt would be as the unit potential flow streamline vector on 
the body. However, streamlines on a body in bounded and unbounded flows, 
described by © and oR; respectively, do not coincide in general, except in 
the special case of an axisymmetric body in a flow facility of circular 
cross section. In the case of a ship hull, if T= (1,0,0), and 55 = 

2° L * T under the assumption that the ship is thin, Equation (4) can be 


reduced to 


Au=u-u (5) 


where 


0 
2c ar |¢ (3 -¥2) - @(-F>¥°) Je 
-T 
0 
uy “An : T { E (4 .y,0) oC (- ¥ .y,0] oy 
-T 


where L and T are the ship length and draft, respectively. In Equation (5), 
the draft T is assumed to be uniform from the bow at x =- L/2 to the stern 
at x = L/2; the centerplane of the ship is on z = 0. 

Similarly, for a slender axisymmetric body of revolution in a wind 
tunnel of circular cross section, the mean-speed increment averaged over 


the body surface is given by 


fe) 
where 
L — 
eT bei oy =0 
hell) L 
See ISTE R=0 
(6) 
ta 1 x= a =0 
Ws saa 
3 E x =- 3 > R=0 
where 
2 


and the peripheral length along a body meridian is approximated by the 


body length, assuming that the body is slender. 


APPROXIMATE MEAN-SPEED INCREMENT 

In this subsection, the method of obtaining an approximate speed 
correction formula is given, based on the potential jump discussed earlier 
by Ballo Define K as a jump in the velocity potential $ given in 
Equation (1) between the infinite upstream and downstream directions. The 
potential jump K is given by integrating the speed increment along a line 
in the fluid from a point infinitely far upstream to a point infinitely 
far downstream. Numerical solutions for practical ship forms at sub- 
critical speeds in towing tanks and for slender bodies of revolution in 
wind tunnels indicate that most of the potential jump occurs along the 
body length. This finding, observed in numerical solutions, will be used 
as the basis for obtaining the present approximate formula for the speed 
correction. It is possible to prove this empirical finding by showing 
that the values of the potential at the upstream and downstream stagnation 
points are approximately equal to the corresponding asymptotic values of 
the potential in the simple case of axisymmetric flow. However, the proof 
will not be discussed here. Thus, the mean speed increment u due to 


blockage is approximated by 


(7) 


Cc 
G 
i] 
lA 


In a recent simple analysis,* the expressions for the potential jump 
K in terms of the effective volume and the depth Froude number Fy in three 


dimensions with a free surface were 


_ (# +m'/o) U 


2. 
WH (1-F,, ) 


K (8) 


*A more detailed analysis in general cases has been submitted in a 
paper to the Journal of Fluid Mechanics (1978). 


where ¥ = displaced volume 
0 = density of water 
= towing speed 


U 

W = tank depth 
H = water depth 
m 


= m' (F.,) = added mass in the longitudinal direction 


F., = U/vgH 


In the derivation of Equation (8), it is assumed that the waterplane area 
of the ship hull is so thin that a line integral term is neglected, and the 
body boundary condition is satisfied exactly on the body surface. From 


Equations (7) and (8) is obtained 


Asi) SP cel y/o 


U 2 
NG KGa) 


(9) 


where A = WH is the cross sectional area of the tank. It is of interest 
to note that when the value of g approaches infinity, Fu approaches zero, 
and Equation (9) reduces to the case of a wind tunnel, where 

Au id Mest: m'/0o 


‘um AL (10) 


It is also of interest to note that when the body boundary condition is 
linearized, i.e., satisfied on the body centerplane, Equations (9) and (10) 


further reduce to 


Ta eee eo 
AL (1-F, ) 


in the presence of a free surface and 


Aaa a ' 
UATE ao") 
in the absence of a free surface. 
APPLICATIONS 


TOWING TANK EXPERIMENT 

To test the new blockage correction formula, three sets of compu- 
tations were first made for the same model in three different towing tanks. 
The first two tanks had the dimensions given by Tamercaea es see Table l. 
The third tank was approximately four times greater in cross sectional area 
than the large tank listed in Table 1, i.e., W = 24 m and H=12m. The 
specific ship model considered was the Wigley parabolic model (Model M1719 


in Tamura), and the equation of the hull surface was given by 


2 


N 
i] 
I+ 
hol| 


2 
tee | 1-(2) (11) 
where L/B = 10, and T/L = 0.0625. The geometric particulars of the models 
have been given in Table 2. 

In the computations, the ship hull boundary condition was linearized; 
thus, speed correction formula (Equation (9')) was used. To test the 
present mean-speed correction formula, computations were also made from 
Equation (5) the exact mean-speed increment averaged over the hull surface 
from the local velocities obtained by the finite-element methods In 
computing the value of u, from Equation (5), the numerical result for the 
extra large tank was used in place of the perturbation potential for un- 
bounded water Oe because the effect of the tank wall and the bottom was 
found to be negligibly small. Comparisons between the "exact" and 
approximate mean-speed increments are given in Table 3. Agreement is 
reasonably good. It should be noted in Table 3 that the exact mean speed 
averaged on the hull surface us defined by Equation (5), is not only 


nonzero but also independent of Froude number. It should also be noted 


TABLE 1 - DIMENSIONS OF SMALL, LARGE, AND EXTRA LARGE 


Extra Large 
Tank 
24 


TOWING TANKS"? 2 


Width in meters 
Mean Water Depth in meters 182) 


TABLE 2 - WIGLEY PARABOLIC MODEL (TAMURA MODEL M1719) 


Length between Perpendiculars in meters 8.000 
Length of Waterline in meters 7.984 
Beam in meters 0.800 
Draft in meters 0.500 
Volume in meters 1.422 
Wetted Surface in meters 9.408 
Block Coefficient 0.4453 


Prismatic Coefficient 0.6680 


TABLE 3 — COMPARISONS OF MEAN-SPEED INCREMENT, COMPUTED BY 
NUMERICAL RESULTS AND BY PRESENT FORMULA FOR 


= * 
FS 0.4 


Exact Numerical Results Equation (9') 
Tank a J i Lad a 
u/U u/U Au/U = (u-u)/U Au/U 


Small} 0.017198 | 0.030514 0.0133 0.0128 
Large | 0.017198 | 0.019425 0.0022 0.0029 


*Results of extra large tank were used to compute Uo 
as discussed in text. 


that the free surface effect on the velocity profile on the body surface 
would be significantly dependent upon whether the hull is in a shallow 

towing tank or in unbounded water. The present study indicates that the 
approximate speed correction formula satisfactorily treats the seemingly 
complicated free-surface effect on the mean-speed increment on the body. 


The total resistance coefficient C determined experimentally by 


EY 
Tamura, and the wave resistance coefficients Cx computed by the finite 
element method, are given in Table 4. 

In presenting our results, the total resistance coefficient Co and 


the wave resistance coefficient Ci are defined as 


2s Dee 2/3 
C= /e 


(12) 
ch ai eee 
Ww w/ 2 
where Ro and R are, respectively, the total and wave resistances. 
The frictional resistance coefficients, Cy and cu are defined by 
CP as u SS 
S (13) 
AGE Pa peaiobny 2 (Ole e2 22/3 
eae URE sais 7/38 ¥ 


where Sy is the model wetted surface area. The model length Froude number 


F. and Reynolds number RO are defined by 


I; 
FL = U/¥gL 
(14) 
Re v/VOL 


10 


where V is the kinematic viscosity of water. Here the Reynolds number R 
is obtained by assuming that the freshwater temperature in two towing 
tanks was 20 C. Table 4 results are given in Figure 1; the wave resistance 
computed for the extra large tank was taken to be the same as for unbounded 
water, already mentioned. In Figure 1, hull wave resistance in the large 
tank is very close to that for the extra large tank. Thus, the blockage 
effect on wave resistance is very small for the large tank. Also, the 
main difference in the total resistance coefficients Cy measured in the 
small and large towing tanks is due primarily to the difference in the 
model wave resistance computed for the two tanks. 

Table 4 gives the speed corrections computed from Equation (9') along 
with the corrected values of (C,-C_). The corrected value of (C.-C) is 
given by (C,-C ) (u/ (uU+Au)). Table 5 gives the frictional resistance co- 
efficients Ch and Co» computed from International Towing Tank Conference 
(ITTC) (1957) and American Towing Tank Conference (ATTC) friction formulas. 
In the present study, it is assumed that the total resistance less the 
computed theoretical wave resistance is approximately equal to the 
frictional resistance, since the ship hull is thin and smooth, i.e., form 
drag is assumed to be negligibly small. If we make use of the @renailile 
correlation of partial form factor k. with hull-fineness parameter 
= C. ¥(B/L) (2T/L) for the Wigley parabolic model with Fh = 0.5, we find 
that k = 0.04; i.e., form drag is estimated to be only 4 percent of the 
frictional drag and a still lower percentage of the total drag. Accord- 
ingly, speed correction Equation (9') was applied to the resistance 
component co = Co = ce to correct for blockage effect. Results given in 
Table 5 are shown in Figure 2. In Figure 2, the corrected values of 
(C.-C) are lower than the values of CE given by ITTC and ATTC friction 
formulations, indicating negative hull form drag, which is not acceptable. 
In other words, if the form drag coefficient and other corrections had 
been added to the values of ITTC and ATTC friction coefficients, this 


discrepancy would be even larger. The discrepancy seems to have been 


caused by computed values of the wave resistance being too large. 


a 


LEGEND: 


SMALL TANK 
———=; LARGE TANK 
re) EXTRA LARGE TANK 


RESISTANCE COEFFICIENTS 


0.325 0.35 0.375 0.40 


F, = U/JoL 


Figure 1 - Resistance Coefficients, Cos Cis and 
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C 
F 


TABLE 4 - RESISTANCE COEFFICIENTS OF THE WIGLEY PARABOLIC MODEL AT 
TWO DIFFERENT TOWING TANKS (MODEL M1719) 


Cc C 
a Ww (C.-C) (FL) 


Experiment Numerical 
by Tamura by Bai 


Corrected Corrected 


SO OO SHO) CO © © 


Values in both tanks are converted to the case without blockage. 


TABLE 5 - FRICTIONAL RESISTANCE COEFFICIENTS Ch AND cm 


COMPUTED FROM ITTC (1957) AND ATTC FRICTION FORMULAS, 
AT A FRESHWATER TEMPERATURE OF 20 C 


DEAN; ALE) 3y7/ ATTC 
F R 
L n A A 


0.003565 0.003444 
0.003516 0.003400 
0.003470 0.003360 
0.003429 0.003323 
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A 
ITTC 1957 Ce 


° Ee ° 


wees . 


0.325 0.35 0.375 0.40 0.425 
FL =U VoL 


Figure 2 - Corrected Values of (Cp-C_) from Small Tank 
o, Large Tank x, and ITTC 1957 and ATTC Curves C. 
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In the numerical computation of wave resistance by the finite-element 
method, 44 nodes on the ship hull surface, i.e., on the centerplane, and 
1496 nodes for the entire fluid domain were taken. One may expect more 
refined results by reducing the size of finite elements. To treat low 


values of Froude number accurately, smaller and more elements are necessary. 


WIND TUNNEL EXPERIMENT 

As a second example, the blockage effect was considered for a wind 
tunnel having a uniform circular cross section of radius Ro: The specific 
body geometry considered was a prolate spheroid with its meridian profile 


given by 


silt (15) 


for the special case when a/b = 4. 

The potential flow for the axisymmetric boundary configurations con- 
sidered herein could have been computed by the conventional method of 
integral equations; i.e., the axial source and doublet distributions or 
the vortex sheet on the surface, etc., as discussed in Landweber. 

However, the velocity potential has been computed by the finite-element 
method. Computations have been made for seven values of Ro/b salen 25) 

655 25 Sy45 54 nach IS) eulilsatopev y/o a, \eiovexal R/b = 15 was computed, the 
effect of the tunnel wall on the body surface was negligibly small as if 
the body were moving in an infinite fluid. The value of u,/U defined in 
Equation (5), computed by using the result of Ri/b = 15, was 0.08185, 
whereas that computed by using the exact analytic result for the unbounded 
Water, -Laec, Ri/b =o, given in item was 0.08156. 

The computed velocity potential > is shown in Figure 3 for three 
values of Ri/b = 1.25, 1.5, and 15. To illuminate the assumption made to 
obtain the present approximate mean-speed correction, Figure 3 shows 


straight lines drawn from the origin to the asymptotic values of K/2 at the 


15 


¢/Ub 


LEGEND: 


ON BODY 

— — — AT R/b = 1.25 

—.—— LINEAR POTENTIAL VARIATION 
ASSUMED IN THE PRESENT SPEED 
CORRECTION FORMULA (THE 
SLOPE IS THE SPEED CORRECTION) 


Figure 3 - Velocity Potential for a Spheroid (a/b=4) 
in a Wind Tunnel with a Circular Cross Section 
of Radius Ro 
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the downstream stagnation point x = L/2. The slope of each straight line 
is equal to the speed correction defined by Equation (7). Owing to the 
skew symmetry of the potential with respect to x = 0, the result for the 
upstream half-body can be obtained from the downstream potential shown in 
Figure 3. The velocity potential increases monotonically from a value 
slightly lower than - K/2 at the upstream stagnation point to a value 
slightly higher than K/2 at the downstream stagnation point on the body 
surface. However, the potentials at R = 1.25b approach monotonically the 
asymptotic values at both ends for Ri/b S 1525 aml R,/b St 355 

In Table 6 the approximate mean speed correction given by Equation 
(10) is compared with the exact mean-speed correction computed from 
Equation (6). Table 6 also gives the speed correction obtained by the 


Lock and Johansen formula, which is given in Ponee as 
= 3 
ta = 2,501 iE ) (16) 


When Ri/b < 3, our approximate results show better agreement with the exact 
numerical results than with those of Lock and Johansen. 

In Figure 4 computed values of the added mass coefficient and the 
mean speed correction Au/U are shown as a function of b/R.- In Figure 4, 
note that for b/R, > 0.765, the contribution of the added mass to the 
speed correction in Equation (10) is more dominant than the contribution 
of the displaced volume, i.e., m=m'/p¥ > 1. This finding indicates that 
a crude blockage correction, based on only the local cross sectional area 
of the body using one-dimensional theory, cannot always give a good 
approximation of the mean-speed correction when the added mass coefficient 


is not small. 
CONCLUSIONS 


In the present study a new mean-speed formula for corrections caused 


by blockage is proposed. The approximate formula is tested by comparing 


Ly) 


m AND Au/U 


TABLE 6 — COMPARISONS OF MEAN-SPEED INCREMENTS ON A 
SPHEROID IN A WIND TUNNEL, COMPUTED BY A NUMERICAL 
METHOD, APPROXIMATE FORMULAS (u ,/U=0. 0813557 


OBTAINED BY LAMB WAS USED) 


Present 
Formula 
(Equation (6)) 


Lock and Johansen 
(Equation (16)) 


1.25 


i aes | 


b/R,, 


Figure 4 - Added Mass Coefficient m and Speed Correction 
Au/U for a Spheroid in a Circular Wind Tunnel 
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it with an exact numerical mean-speed correction, computed by the finite- 
element method for both a towing tank experiment and a wind tunnel experi- 
ment. The two predictions are shown to be in good agreement for both 
facilities. It is shown that the effect of added mass coefficient on the 
speed correction of a body is very significant as the blockage effect in- 
creases. It is also found that the main difference in the total resistance 
coefficient measured in a large and a small towing tank is due primarily 

to the difference in the model wave resistances computed for the two tanks. 
Further investigation is necessary to take into account other blockage 
corrections due to viscous effects such as flow separation and wake dis- 


placement thickness effects. 
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